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To date, detailed core-collapse supernova simulations have not succeeded in elu- 
cidating the explosion mechanism. But whatever the mechanism turns out to be, 
there will be a neutrino- heated outflow between the hot, newly-born neutron star 
and the outgoing supernova shock wave whose neutron richness, low density, and 
high temperature could provide promising conditions for the synthesis of heavy 
nuclei via rapid neutron capture (the r-process). General relativistic effects im- 
prove the prospects for this outflow as an r-process site. Among relativistic effects, 
the enhanced "gravitational potential" is more important than the gravitational 
redshift or trajectory bending of neutrinos. 



1. Core-collapse Supernovae 

Core-collapse supernovae — those of Type lb, Ic, and II — result from the 
catastrophic collapse of the core of a massive star. For most of their exis- 
tence, stars burn hydrogen into helium. In stars at least eight to ten times 
as massive as the Sun, temperatures and densities become sufficiently high 
to burn to carbon, oxygen, neon, magnesium, and silicon and iron group 
elements. The iron group nuclei are the most tightly bound, and here burn- 
ing in the core ceases. The iron core — supported by electron degeneracy 
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pressure — eventually becomes unstable. Its inner portion undergoes ho- 
mologous collapse (velocity proportional to radius), and the outer portion 
collapses supersonically. Electron capture on nuclei is one instability lead- 
ing to collapse, and this process continues throughout collapse, producing 
neutrinos. These neutrinos escape freely until densities and temperatures 
in the collapsing core become so high that even neutrinos are trapped. 

Collapse is halted soon after the matter exceeds nuclear density; at this 
point (called "bounce"), a shock wave forms at the boundary between the 
homologous and supersonically collapsing regions. The shock begins to 
move out, but after the shock passes some distance beyond the surface of 
the newly born neutron star, it stalls as energy is lost to neutrino emission 
and dissociation of heavy nuclei falling through the shock. 

The details of how the stalled shock is revived sufficiently to continue 
plowing through the outer layers of the progenitor star are unclear. Some 
combination of neutrino heating of material behind the shock, convection, 
instability of the spherical accretion shock, rotation, and magnetic fields 
launches the explosion. 

It is natural to consider neutrino heating as a mechanism for shock 
revival, because neutrinos dominate the energetics of the post-bounce evo- 
lution. Initially, the nascent neutron star is a hot thermal bath of dense 
nuclear matter, electron/positron pairs, photons, and neutrinos, contain- 
ing most of the gravitational potential energy released during core collapse. 
Neutrinos, having the weakest interactions, are the most efficient means of 
cooling; they diffuse outward on a time scale of seconds, and eventually 
escape with about 99% of the released gravitational energy. 

Because neutrinos dominate the energetics of the system, a detailed un- 
derstanding of their evolution will be integral to any detailed and definitive 
account of the supernova process. If we want to understand the explosion — 
which accounts for only about 1% of the energy budget of the system — we 
should carefully account for the neutrinos' much larger contribution to the 
energy budget. 

What sort of computation is needed to follow the neutrinos' evolution? 
Deep inside the newly-born neutron star, the neutrinos and the fluid are 
tightly coupled (nearly in equilibrium) ; but as the neutrinos are transported 
from inside the neutron star, they go from a nearly isotropic diffusive regime 
to a strongly forward-peaked free-streaming region. Heating of material be- 
hind the shock occurs precisely in this transition region, and modeling this 
process accurately requires tracking both the energy and angle dependence 
of the neutrino distribution functions at every point in space. 
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A full treatment of this six-dimensional neutrino radiation hydrody- 
namics problem remains too costly for currently available computational 
resources. Throughout the 1990s, several groups performed simulations in 
two spatial dimensions with simplified neutrino transport. 

One simplification allowed for neutrino transport in two spatial di- 
mensions, but with neutrino energy and angle dependence integrated 
out — effectively reducing a five dimensional problem to a two dimensional 
one. 1 ' 2 ' 3 These simulations showed convection in two regions. First, loss of 
electron neutrinos from the outer layers of the neutron star caused composi- 
tion gradients that could drive convection, which boosted neutrino luminosi- 
ties by bringing hotter material to the surface. Second, heating decreased 
further from the neutron star surface, giving rise to a negative entropy gra- 
dient. The resulting convection increased the efficiency of neutrino heating 
by delivering heated material to the region just behind the shock. These 
simulations exhibited explosions, suggesting that the enhancements in neu- 
trino heating behind the shock resulting from convection provided a robust 
explosion mechanism. More recent simulations in three spatial dimensions 
with this same approximate treatment of neutrino transport showed similar 
outcomes. 4 

A different simplification of neutrino transport employed in the f990s 
was the imposition of energy-dependent neutrino distributions from spher- 
ically symmetric simulations onto fluid dynamics computations in two spa- 
tial dimensions. 5 Unlike the multidimensional simulations discussed above, 
these did not exhibit explosions, casting doubt upon claims that convection- 
aided neutrino heating constituted a robust explosion mechanism. 

The nagging qualitative difference between multidimensional simula- 
tions with different neutrino transport approximations renewed the moti- 
vation for simulations in which both the energy and angle dependence of the 
neutrino distributions were retained. Of necessity, the first such simulations 
were performed in spherical symmetry (actually a three-dimensional prob- 
lem, depending on one space and two momentum space variables) . Results 
from three different groups are in accord: Spherically symmetric models do 
not explode, even with solid neutrino transport. 6 ' 7 ' 8 

Recently, one of these groups performed simulations in two spatial di- 
mensions, in which their energy- and angle-dependent neutrino transport 
was made partially dependent on spatial polar angle as well as radius. 9,10 
Explosions were not seen in any of these simulations, except for one in 
which certain terms in the neutrino transport equation corresponding to 
Dopplcr shifts and angular aberration due to fluid motion were dropped. 
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This was a surprising qualitative difference induced by terms contributing 
what are typically thought of as small corrections. The continuing lesson is 
that getting the details of the neutrino transport right makes a difference. 

Where, then, do simulations aiming at the explosion mechanism stand? 
The above history suggests that elucidation of the mechanism will require 
simulations that feature truly spatially multidimensional neutrino trans- 
port. Development of the formalism, 11 algorithms, 12 and computer code 
necessary to this transport capability is ongoing as part of the Terascale 
Supernova Initiative. This substantial collaboration — led by Anthony Mez- 
zacappa of Oak Ridge National Laboratory and funded by the Department 
of Energy's Scientific Discovery through Advanced Computing (SciDAC) 
program — is dedicated to the elucidation of the core-collapse supernova ex- 
plosion mechanism through supercomputer simulations. At least one other 
major development of multidimensional neutrino transport capability is un- 
derway as well. 13 In addition to sophisticated neutrino radiative transfer, 
inclusion of magnetic field dynamics seems increasingly strongly motivated 
as a possible driver of the explosion, because simulations with "better" 
neutrino transport have failed to explode — even in multiple spatial dimen- 
sions. Work on the inclusion of magnetic fields is also part of the Terascale 
Supernova Initiative. 

2. Relativistic Neutrino-driven Winds: Early Work 

It was realized in the early 1990s that whatever the explosion mechanism 
turns out to be, there will be an evacuated region between the hot, newly- 
born neutron star and the outgoing supernova shock wave, whose neutron 
richness, low density, and high temperature of the region could provide 
promising conditions for the synthesis of heavy nuclei via rapid neutron 
capture (the r-process). 14 Consideration of the exploding supernova mod- 
els a of Wilson and Mayle in the late 1980s and early 1990s 18 led Meyer et 
al. 14 to this insight, and to the related observation that the ejected amount 
of this "hot bubble" material — together with the Galactic supernova rate — 
seemed about right to account for the amount of r-process material in the 



a Thcsc simulations were spherically symmetric. In light of the discussion in the previous 
section, it may be surprising that explosions were obtained. However, a prescription 
mocking up the multidimensional effects of a doubly- diffusive fluid instability (the so- 
called "neutron fingers") was included in these simulations, which boosted neutrino 
luminosities sufficiently for the neutrino-driven explosion mechanism to succeed. That 
the necessary conditions actually exist for this particular instability to operate effectively 
has been disputed, 15,16 ' 17 but related instabilities may produce similar results. 17 
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Galaxy. 

This proposal was buttressed by subsequent work of Woosley et al.: 19 
The Wilson and Mayle models were run out to several seconds past core 
bounce and explosion, showing that the intense neutrino fluxes emitted by 
the cooling neutron star did indeed drive neutron-rich matter off its surface 
into the evacuated region below the shock wave. This work also included 
nucleosynthesis calculations — performed by post-processing matter trajec- 
tories obtained in these simulations — that yielded impressive agreement 
with the observed solar system r-process abundance distribution. 

The observed r-process abundances require that ~100 neutrons be avail- 
able for capture on each iron peak "seed" nucleus, and this neutron/seed 
ratio is determined by three parameters characterizing the astrophysical 
environment: the entropy per baryon S, the electron fraction Y e , and the 
dynamic expansion time scale Td yn - High entropy favors the relative disor- 
der of free nucleons, as opposed to their being locked up in heavier nuclei 
(high rates of photodisintegration at the high temperatures and low densi- 
ties associated with large S provide the microscopic mechanism). Low elec- 
tron fraction corresponds to neutron richness: Y e = n p j (n n + n p ), where n p 
and n n are the number densities of protons and neutrons (including those 
locked up in nuclei). A short dynamic expansion time scale prevents too 
many seed nuclei from building up, by causing the freeze-out of bottleneck 
three-body reactions producing 12 C from 4 He. 

The "hot bubble" exhibited in Woosley et al. 19 appears to be a sub- 
sonic outflow that bumps up against the shock wave sitting about 10 km 
from the neutron star. In these conditions, the factor most favorable to the 
r— process turns out to be a high entropy of S <~ 400 (in units of Boltz- 
mann's constant); values of Y e ~ 0.4 and Td yn ~ls are modest. 

Several factors motivate modeling this neutron-rich outflow from the 
neutron star surface on its own — separate from large-scale supernova 
simulations — in order to gauge its suitability as an r-process site. One 
obvious roadblock to the use of the large-scale simulations is the failure of 
most models to explode, as described in the previous section. Moreover, 
even if explosions are obtained, many of these simulation codes are not 
well suited to the task of running to the late times (> 10 s) needed to 
follow the wind. A physical justification for simple wind models is that ex- 
ploding two-dimensional supernova models with convection settled down to 
approximately spherically symmetric and stationary conditions as the wind 
phase approached. 2,3 Simple models can provide physical insight, and are 
more amenable to parameter studies (e.g. dependence on neutron star mass 
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Figure 1. The final entropy per baryon in units of Boltzmann's constant, as a function 
of the supernova core Schwarzschild radius divided by the core radius. The circle is 
from the Qian and Woosley numerical calculation of model 10B with post-Newtonian 
corrections. 

and radius, neutrino luminosities and average energies) than the large-scale 
simulations. 

Studies of simpler models of the neutrino-driven wind began with an 
important paper of Qian and Woosley. 20 They obtained estimates of S, Td yn , 
and mass outflow rate from analytic "wind" models based on steady-state 
Newtonian fluid equations describing the matter outflow. They showed that 
the putative high entropy 19 of the "hot bubble" was difficult to explain, 
even with an outer boundary pressure. Their values of S fell well short of 
that required for a robust r-process. They confirmed their analytic results 
with a hydrodynamic code that included simple input neutrino heating. In 
a few of their numerical runs, they employed an enhanced "gravitational 
force" — ^*r- — > — 1 _2GM/r mo ti va ted by the corresponding term in the 
relativistic fluid equations in Schwarzschild geometry. These cases resulted 
in larger values of S and smaller values of r^ym an d it was pointed out that 
both went in the right direction towards more favorable conditions for the 
r-process. 

Cardall and Fuller 21 pursued this hint on the effects of relativity, follow- 
ing the Qian and Woosley 20 approach, but with relativistic fluid equations. 
In addition to the relativistic effects in the fluid equations (e.g. enhanced 
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Figure 2. The dynamic expansion time scale as a function of the supernova core 
Schwarzschild radius divided by the core radius. The circle is from the Qian and Woosley 
numerical calculation of model 10B with post-Newtonian corrections. 
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Figure 3. The mass outflow rate as a function of the supernova core Schwarzschild 
radius divided by the core radius. The circle is from the Qian and Woosley numerical 
calculation of model 10B with post-Newtonian corrections. 

"gravitational force," Lorentz factors limiting velocities to the speed of 
light, internal energy density and pressure contributing to inertia), rela- 
tivistic effects on the neutrino heating were taken into account, specifically 
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gravitational redshift and the bending of neutrino trajectories. 

These calculations affirmed that general relativistic effects make con- 
ditions in the wind more hospitable to the r-process. Figures 1, 2, and 3 
show S, Tdyn, and the mass outflow rate as a function of the compactness 
of the neutron star (Schwarzschild radius divided by radius). b Trends of 
increasing S and decreasing Td yn with increasing compactness are both fa- 
vorable, though this comes at a cost of a smaller mass outflow rate, which 
translates into the production of less r-process material. These figures are 
not physically meaningful for 2< ^Jf > 0.66, as compactness greater than 
this would imply an equation of state that violates causality (sound speed 
greater than the speed of light). Comparing Figures 1, 2, and 3 with a pa- 
rameter study 22 of the neutron/seed ratio as a function of S, Td yn , and Y" e , 
Cardall and Fuller 21 concluded that suitable conditions for the r-process 
could just be achieved near the neutron star causality limit, but that the 
amount of ejected material becomes uncomfortably (though perhaps not 
prohibitively) small. 

The most important relativistic effect appears to be the enhanced "grav- 
itational force." One way to see this is to note the good agreement be- 
tween the Cardall and Fuller 21 results and the single displayed Qian and 
Woosley 20 numerical result, in which the enhanced "gravitational force" 
was the only relativistic effect included. This agreement is particularly 
striking in the case of the entropy per baryon S. Roughly speaking, the 
entropy per baryon is the energy a baryon acquires by neutrino heating, 
divided by a characteristic temperature at which the heating takes place. 
If a baryon is to escape the "gravitational potential," the acquired energy 
per baryon is roughly the baryon mass times the gravitational potential at 
the neutron star surface. 20 In the relativistic case, the "gravitational poten- 
tial" obtained from the enhanced "gravitational force" expressed above is 
\ In (l — ^r-), which reduces to the Newtonian expression for small 
compactness. Gravitational neutrino redshift and trajectory bending have 
some effect on the characteristic temperature at which the heating takes 
place, 21 but this effect is modest in comparison with the logarithmic de- 
pendence of the "gravitational potential." The result is that Figure 1 for 
S tracks the logarithmic dependence of the gravitational potential quite 
closely; the corresponding plot for the Newtonian case would be a straight 



This variation was actually computed by varying the neutron radius with the mass held 
fixed at IAMq. A particular prescription for the variation of neutrino luminosity was 
also employed. 
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line matching the relativistic curve at low compactness. It is evident that 
relativity makes a nontrivial difference for neutron star masses and radii. 

3. Recent Work and Outstanding Issues 

Cardall and Fuller 21 concluded that the prospects for suitable r-process 
conditions in the neutrino-driven wind improve from something like 'rather 
pessimistic' in the Newtonian case 20 to perhaps 'not inconceivable' when 
general relativity is taken into account — an assessment confirmed in subse- 
quent work. 23 ' 24 ' 25,26 These works involved various improvements on pre- 
vious semi-analytic estimates, including full numerical solution of the wind 
equations; comprehensive variations of neutron star mass, radius, and neu- 
trino luminosities; tracking of the electron fraction; construction of evo- 
lutionary sequences; proper treatment of transonic winds; and r-process 
network nucleosynthesis calculations. The consensus emerging from these 
works was that suitable r-process conditions might obtain, but that it would 
be in the form of a modest entropy, rapidly expanding (and possibly tran- 
sonic) "wind" rather than a high entropy, subsonic "hot bubble." The 
rather massive and compact neutron stars that seem to be required do not 
seem likely given current understanding of the dense nuclear matter equa- 
tion of state, but new analyses of neutrino interactions with this dense mat- 
ter may alleviate this problem by implying higher neutrino luminosities. 27 

A verdict of 'not inconceivable' is not particularly comforting, but given 
the apparent preference of Galactic chemical evolution models for a super- 
nova r-process source over neutron star binary mergers, 28 fresh ideas on 
supernova winds are welcome. Hydrodynamic calculations have been per- 
formed in which outer boundary effects are claimed to play a significant 
role, allowing suitable r-process conditions for neutron stars of canonical 
mass; 29 perhaps the issue of "rapid wind" vs. "hot bubble" is worth an- 
other examination. A very different and interesting idea is that magnetic 
fields may trap matter long enough to be significantly heated before being 
released into the wind. 30 ' 31 

Final understanding of the suitability of the supernova neutrino-driven 
wind as an r-process site will probably only come in the context of an under- 
standing of the explosion mechanism — which brings us back to the discus- 
sion of large-scale simulations in the first section. Once this understanding 
is achieved through detailed simulation, we will have better qualitative and 
quantitative understandings of neutrino heating, magnetic fields, fallback 
of material at late times, and other phenomena that will influence the wind. 
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